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Abstract — In this work it is shown how to obtain, in a simple way, localized (non- 
diffractive) subluminal pulses as exact analytic solutions to the wave equations. These 
q new ideal subluminal solutions, which propagate without distortion in any homogeneous 

linear media, are herein obtained for arbitrarily chosen frequencies and bandwidths, 
avoiding in particular any recourse to the non-causal components so frequently plaguing 
the previously known localized waves. The new solutions are suitable superpositions of 
— zeroth-order, in general — Bessel beams, which can be performed either by integrating 
with respect to (w.r.t.) the angular frequency u, or by integrating w.r.t. the longitudinal 
wavenumber k z : Both methods are expounded in this paper. The first one appears 
to be powerful enough; we study the second method as well, however, since it allows 
dealing even with the limiting case of zero-speed solutions (and furnishes a new way, in 
terms of continuous spectra, for obtaining the so-called "Frozen Waves", so promising 
also from the point of view of applications). We briefly treat the case, moreover, of non 
axially-symmetric solutions, in terms of higher order Bessel beams. At last, particular 
attention is paid to the role of Special Relativity, and to the fact that the localized waves 
are expected to be transformed one into the other by suitable Lorentz Transformations. 
The analogous pulses with intrinsic finite energy, or merely truncated, will be constructed 
in another paper. In this work we fix our attention especially on electromagnetism and 
optics: but results of the present kind are valid whenever an essential role is played by a 
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wave-equation (like in acoustics, seismology, geophysics, gravitation, elementary particle 
physics, etc.). 
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1 Introduction 

Since more than ten years, the so-called (non-diffracting) "Localized Waves" (LW), which 
are new solutions to the wave equations (scalar, vectorial, spinorial,...), are in fashion, 
both in theory and in experiment. In particular, rather well-known are the ones with 
luminal or superluminal peak- velocity: Like the so-called X-shaped waves (see [221 [38] an d 
refs. therein; for a review, see, e.g., ref. [43]), which are supersonic in Acoustics |24j, and 
superluminal in Electromagnetism (see [23] and refs. therein). 

Since Bateman[20] and later on Courant & Hilbert[T9], it was already known that 
luminal LWs exist, which are as well solutions to the wave equations. More recently, 
some attention[31 ESI EDI Ell E] started to be paid to subluminal solutions too. Let 
us recall that all the LWs propagate without distortion — and in a self-reconstructive 
way [551 121 EH] — in a homogeneous linear medium (apart from local variations). 

Like in the superluminal case, the (more orthodox, in a sense) subluminal LWs can 
be obtained by suitable superpositions of Bessel beams. They have been till now almost 
neglected, however, for the mathematical difficulties met in getting analytic expressions 
for them, difficulties associated with the fact that the superposition integral runs over a 
finite interval. We want here to re-address the question of such subluminal LWs, showing, 
by contrast, that one can indeed arrive at exact (analytic) solutions, both in the case of 
integration over the Bessel beams' angular frequency u, and in the case of integration over 
their longitudinal wavenumber k z . The first approach, herein investigated in detail, is 
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enough to get the majority of the desired results; we study also the second one, however, 
since it allows treating the limiting case of zero-speed solutions (and furnishes a new 
method — based on a continuous spectrum — for obtaining the so-called "Frozen Waves" , 
so promising also from the point of view of applications). Moreover, we shall deal with non 
axially-symmetric solutions, in terms of higher order Bessel beams. At last, particular 
attention is paid to the role of Special Relativity, and to the fact that the localized waves 
are transformed one into the other by Lorentz Transformations. 

The present article is devoted to the exact, analytic solutions: i.e., to ideal solutions. In 
another paper, we shall go on to the corresponding pulses with finite energy, or truncated, 
sometimes having recourse — for simplicity, and only in such more realistic cases — to 
approximations. We shall fix our attention especially on electromagnetism and optics: 
but our results are valid whenever an essential role is played by a wave-equation (like in 
acoustics, seismology, geophysics, gravitation, elementary particle physics, etc.). 

Let us recall that, in the past, too much attention had not been paid even to 1983 
Brittingham's paper [41], where he claimed the possibility of obtaining pulse-type solu- 
tions to the Maxwell equations, which propagate in free space as a new kind of speed-c 
"solitons" . This was partially due to the fact that that author was unable to get correct 
finite-energy expressions for such wavelets, and to make suggestions about their practi- 
cal production. Two years later, however, Sezginerjl2] was able to show how to obtain 
quasi-nondiffracting luminal pulses endowed with a finite energy. Finite-energy pulses do 
no longer travel undistorted for an infinite distance, but they can nevertheless propagate 
without deformation for a long field-depth, much larger than the one achieved by ordinary 
pulses, as the gaussian ones: Cf., e.g., refs. [25J E3 EHl EH El E] and refs. therein. 

Only after 1985 the general theory of all LWs started to be extensively developed [3T| 
E21 E31 1221 LI21 12HI EH EH S3 SS] both in the case of beams, and in the case of pulses. 
For reviews, see for instance the refs. [331 HH El [HI |9] and references therein. For the 
propagation of LWs in bounded regions (like guides), see refs. [521 EH E31 ES] and refs 
therein. For the focusing of LWs, see refs .[541 18] and refs therein. As to the construction 
of LWs propagating in dispersive media, cf . refs. [251 SHI SSI E7J [13, [TJ EOj SHI E] ; and, 
for lossy media, see ref.[59j and refs. therein. Al last, for finite-energy, or truncated, 
solutions see refs. [Ml El ESI EU ESI EI] ■ 

By now, the LWs have been experimentally produced [2U EE! EZ], and are being applied 
in fields ranging from ultrasound scanning[28l [291 EH] to optics (for the production, e.g., 
of new type of tweezers [62]). All these works demonstrated that nondiffracting pulses 
can travel with an arbitrary peak-speed v, that is, with < v < oo; while Brittingham e 
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Sezginer had confined themselves to the luminal case (v = c) only. As already commented, 
the superluminal and luminal LWs have been, and are being, intensively studied, whilst 
the subluminal ones have been neglected: Almost all the few papers dealing with them had 
till now recourse to the paraxial [58] approximation [IT], or to numerical simulations |14j. 
due to the above mentioned mathematical difficulty in obtaining exact analytic expressions 
for subluminal pulses. Actually, only one analytic solution is known P, [12j [63J, [301 E] 
biased by the physically unconvenient facts that its frequency spectrum is very large, it 
doesn't even possess a well-defined central frequency, and, even more, that non-causal (i.e., 
backwards traveling) [521 [H] components are a priori needed for constructing it. Aim of 
the present paper is showing how subluminal localized exact solutions can be constructed 
with any spectra, in any frequency bands and for any bandwidths; and, moreover, without 
employing[38l E] any non-causal components. 



2 A first method for constructing physically 
acceptable, subluminal Localized Pulses 

Axially symmetric solutions to the scalar wave equation are known to be superpositions of 
zero-order Bessel beams over the angular frequency u and the longitudinal wavenumber 
k z : i.e., in cylindrical co-ordinates, 

*{p,z,t) = J o duj^dk z S(uj,k z )J \p^—-kl\e^ z e-^ t , (1) 

where k 2 p = u 2 /c 2 — k 2 z is the transverse wavenumber; quantity k 2 has to be positive if 
one does not want to deal with evanescent waves. 

The condition characterizing a nondiffracting wave is the existence [HI [ED] of a linear 
relation between longitudinal wavenumber k z and frequency to for all the Bessel beams 
entering superposition (1); that is to say, the chosen spectrum has to entail [381 EH] for 
each Bessel beam a linear relationship of the typd^ (cf. Figjl]): 

uj = v k z + b (2) 



*More generally, as shown in ref. [35], the chosen spectrum has to call into the play, in the plane oj, k z , 
if not exactly the line (2), at least a region in the proximity of a straight-line of such a type. In the latter 
case, the corresponding solution will possess finite energy, but will have a finite "depth of field" : that is, 
it will be nondiffracting only till a certain finite distance. 
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with b > 0. Requirement (2) can be regarded also as a specific space-time coupling, 
implied by the chosen spectrum S. Equation (2) has to be obeyed by the spectra of any 
one of the three possible types (subluminal, luminal or superluminal) of nondiffr acting 
pulses. Let us mention, incidentally, that with the choice in Eq.(2) the pulse re-gains 
its initial shape after the space-interval Azi = 2nv/b] the more general case can be, 
however, considered [38, 6] when b assumes any values b m = mb (with m an integer), and 
the periodicity space-interval becomes Az m = Azi/m . 

In the subluminal case, of interest here, the only exact solution known till now, 
represented by Eq.(10) below, is the one found by Mackinnon[3j. Indeed, by taking into 
explicit account that the transverse wavenumber k p of each Bessel beam entering Eq.(l) 
has to be real, it can be easily shown (as first noticed by Salo et al. for the analogous 
acoustic solutions [H]) that in the subluminal case b cannot vanish, but must be larger 
than zero: b > 0. Then, by using conditions (2) and b > 0, the subluminal localized 
pulses can be expressed as integrals over the frequency only: 

ty(p,z,t) = exp [— ib-\ / da; SiuS) Jo(pk p ) exp \iuj-] , (3) 

V Jul- V 

where now 

k p = ^2bu-b 2 -(l-v 2 /c 2 )uj 2 (4) 

with 



C = z - vt (5) 

and with 

f b 

1 + v/c 

(6) 

b 

As anticipated, the Bessel beam superposition in the subluminal case results to be an 
integration over a finite interval of u, which, by the way, does clearly shows that the 
non-causal (backwards traveling) components correspond to the interval u_ < ui < b. 
It could be noticed, incidentally, that Eq.(3) does not represent the most general exact 
solution, which on the contrary is a sumJO] of such solutions for the various possible 
values of b just mentioned above: That is, for the values b m = mb and spatial periodicity 
Az m = Azi/m; but we can confine ourselves to solution (3) without any real loss of 
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generality, since the actual problem is evaluating in analytic form the integral entering 
Eq.(3). For any mathematical and physical details, see ref.[B]. 
Now, if one adopts the change of variables 



b /-, v \ 
1 — v z /c z c 



equation (3) becomes 



(7) 



V(p,z,t) = ^-^ exp[-i-z] exp 

c 1 — v z /c z V 



J ds S(s) J 



yjl ~ V 2 /C 2 



v 1 — v 2 /c 2 

b 



exp 



z ~i 27^ S 

c 1 — f z /C z 



(8) 



In the following we shall adhere to some symbols standard in Special Relativity (since the 
whole topic of subluminal, luminal and superluminal LWs is strictly connected [131 123} I2T] 
with the principles and structure of special relativity [cf. [TJ dH] and refs. therein], as we 
shall mention in the Conclusions); namely: 
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7 = 



(9) 



c Vi^W ' 

As already said, Eq.(8) has till now yielded one analytic solution for S(s) = constant, 
only (for instance, S(s) = 1); which means nothing but S(u) = constant: in this case one 
gets indeed the Mackinnon solution [21 [121 E0J 1ST] 



b 2 
2-vy exp 
c 



i- P'J 2 77 

c 



! b 2 



x sine 



7 2 (P 2 + 7 2 C 2 ) 



(10) 



which however, for its above-mentioned drawbacks, is endowed with little physical and 
practical interest. In Eq.(9) the sine function is defined as 



sine x = (sin x)/x 



and 
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c 2 

n = z - Vt, with V = — , (11) 

v 

where V and v are related by the de Broglie relation. [Notice that \& in Eq.(lO), and in 
the following ones, is eventually a function (besides of p) of z, t via ( and 77, both functions 
of z and t ] . 

We can, however, construct by a very simple method new subluminal pulses cor- 
responding to whatever spectrum, and devoid of non-causal (entering) components, by 
taking also advantage of the fact that in our equation (8) the integration interval is finite: 
that it, by transforming it in a good, instead of a harm. Let us first observe that Eq.(lO) 
doesn't admit only the already known analytic solution corresponding to S = constant, 
and more in general to S(cu) = constant, but it will also yield an exact, analytic solution 
for any exponential spectra of the type 

^/ x ,i2nnuj, , . 

S ( U ) = exp [—^-] , (12) 

with n any integer number, which means for any spectra of the type S(s) = 
exp [inn/ (3] exp [irars], as can be easily seen by checking the product of the various expo- 
nentials entering the integrand. In Eq.(12) we have set 



Q = uj + — oj- 
The solution writes in this more general case: 

.b 



#(p,C,77) = 2b(3 Y exp 



i- n 
c 



x exp [ m ~p\ smc ^ 



,; ' 2 7 V+ { h -l 2 C + nix\ . (13) 



c 2 \c 



Let us explicitly notice that also in Eq.(13) quantity n is defined as in Eqs.(ll) above, 
where V and v obey the de Broglie relation vV = c 2 , the subluminal quantity v being 
the velocity of the pulse envelope, and V playing the role (in the envelope's interior) of a 
superluminal phase velocity. 
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The next step, as anticipated, consists just in taking advantage of the finiteness of 
the integration limits for expanding any arbitrary spectra S(u) in a Fourier series in the 
interval w_ < uj < uj + : 

S(u) = A n exp[+in— uj] , (14) 

n=— oo 

where (we can go back, now, from the s to the uj variable): 

' [' u) + 2,71 

A n = — J duS(uj) exp[-m— uj] (15) 

quantity Q being defined as above. 

Then, on remembering the special solution (13), we can infer from expansion (14) that, 
for any arbitrary spectral function S(u), a rather general, axially-symmetric, analytic 
solution for the subluminal case can be written 



V(p,(,v) = 26/3 7 2 exp 



%- (3 T] 

c 



x ^2 A n exp [in—] sine < 
n=— 00 P \ 



^tV+(Vc + J , (16) 



in which the coefficients A n are still given by Eq.(15). Let us repeat that our solution is 
expressed in terms of the particular equation (13), which is a Mackinnon-type solution. 

The present approach presents many advantages. We can easily choose spectra local- 
ized within the prefixed frequency interval (optical waves, microwaves, etc.) and endowed 
with the desired bandwidth. Moreover, as already said, spectra can now be chosen such 
that they have zero value in the region cj_ < uj < b, which is responsible for the non-causal 
components of the subluminal pulse. 

Let us stress that, even when the adopted spectrum S(uj) does not possess a known 
Fourier series (so that the coefficients A n cannot be exactly evaluated via Eq.(15)), one 
can calculate approximately such coefficients without meeting any problem, since our 
general solutions (8) will still be exact solutions. 

Let us set forth in the following some examples. 
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2.1 Some Examples 

In general, optical pulses generated in the lab possess a spectrum centered on some 
frequency value, Uq, called the carrier frequency. The pulses can be, for instance, ultra 
short, when Au/ujq > 1; or quasi monochromatic, when Au/uj << 1, where Auo is the 
spectrum bandwidth. 

These kinds of spectra can be mathematically represented by a gaussian function, or 
functions with similar behavior. 



First Example: 



Let us consider a gaussian spectrum 

S(u) = -^=exp[-a 2 (uj -too) 2 } (17) 
v 7r 

whose values are negligible outside the frequency interval uj_ < u < uo + over which the 
Bessel beams superposition in Eq.(3) is made, it being cj_ = b/(l + /3) and uo + = 6/(1-/3). 
Of course, relation (2) has still to be satisfied, and with b > 0, for getting an ideal 
subluminal localized solution. Notice that, with spectrum (17), the bandwidth (actually, 
the FWHM) results to be Auj = 2/ a. Let us emphasize that, once v and b have been 
fixed, the values of a and u can then be selected in order to kill the backwards-traveling 
components, that exist, as we know (cf. Fig.l) for u < b . 

The Fourier expansion in Eq.(14), which yields, with the above spectral function (17), 
the coefficients 

1 r • 27r i r _n271 " 2 n /-.on 

An ~ w exp [ ~ m li Uo] exp [ ^w^ ] ' ( 8) 



constitutes an excellent representation of the gaussian spectrum (17) in the interval uj- < 
uj < uj + (provided that, as we requested, our gaussian spectrum does get negligible values 
outside the frequency interval < uj < uj + ). 

In other words, on choosing a pulse velocity v < c and a value for the parameter b, a 
subluminal pulse with the above frequency spectrum (17) can be written as Eq.(15), with 
the coefficients A n given by Eq.(18): The evaluation of such coefficients A n being rather 
simple. Let us repeat that, even if the values of the A n are obtained via a (rather good, 
by the way) approximation, we based ourselves on the exact solution Eq.(16). 
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One can, for instance, obtain exact solutions representing subluminal pulses for optical 
frequencies. Let us get the subluminal pulse with velocity v = 0.99 c, angular carrier 
frequency u = 2.4 x 10 15 Hz (that is, A = 0.785 /mi) and bandwidth (FWHM) Au = 
wo/20 = 1.2 x 10 14 Hz, which is an optical pulse of 24 fs (that is the FWHM of the pulse 
intensity). For a complete pulse characterization, one has to choose the value of the 
frequency b: let it be b = 3 x 10 13 Hz; as a consequence one has c<j_ = 1.507 x 10 13 Hz 
and u + = 3 x 10 15 Hz. [This is exactly a case in which the considered pulse is not 
plagued by the presence of non-causal components, since the chosen spectrum forwards 
totally negligible values for u < b}. The construction of the pulse does already result 
satisfactory when considering about 51 terms (—25 < n < 25) in the series entering 
Eq.(16). 

Figures [2] show our pulse, plotted by considering the mentioned fifty-one terms. 
Namely: Fig. (a) depicts the orthogonal projection of the pulse intensity; Fig.(b) shows 
the three-dimensional intensity pattern of the real part of the pulse, which reveals the 
carrier wave oscillations. 




Figure 1: (a) The orthogonal projection of the pulse intensity; (b) The three-dimensional 
intensity pattern of the real part of the pulse, which reveals the carrier wave oscillations. 

Let us stress that the ball-like shap^] for the field intensity should be typically asso- 
ciated with all the subluminal LWs, while the typical superluminal ones are known to be 

'It can be noted that each term of the series in Eq.(16) corresponds to an ellipsoid or, more specifically, 
to a spheroid, for each velocity v. 
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X-shaped[2"2" j |2"5 | |2"T], as predicted, since long, by special relativity in its "non-restricted" 
version: See refs. [U HH1 1231 S3] and refs therein. 



Another interesting spectrum S(u) would be, for example, the "inverted parabola" one, 
centered at the frequency uq: that is, 



S(u) 



-4[^-(^ -A^/2)]^-(^+Aa,/2)] for Uo _ AuJ /2<uj<iJ + ALj/2 

(19) 

otherwise , 



where Aa>, the distance between the two zeros of the parabola, can be regarded as the 
spectrum bandwidth. One can expand S(u), given in Eq.(19), in the Fourier series (14), 
for uo_ < uj < u + , with coefficients A n that — even if straightforwardly calculable — results 
to be complicated, so that we skip reporting them here explicitly. Let us here only mention 
that spectrum (19) may be easily used to get, for instance, an ultrashort (femtasecond) 
optical non-diffracting pulse, with satisfactory results even when considering very few 
terms in expansion (14). 

Second Example: 

As a second example, let us consider the very simple case when — within the integration 
limits c<j_, u> + — the complex exponential spectrum (12) is replaced by the real function 
(still linear in uj) 

S(u) = - ^ rr exp [a(u - cv+] , (20) 

1 — exp [— a(oj + — u)-)\ 

with a a positive number [for a = one goes back to the Mackinnon case] . Spectrum (20) 
is exponentially concentrated in the proximity of u + , where it reaches its maximum value; 
and becomes more and more concentrated (on the left of u + , of course) as the arbitrarily 
chosen value of a increases, with the frequency bandwidth being Auj = 1/a. Let us recall, 
incidentally, that, on its turn, quantity u + and oj- depend on the pulse velocity v and on 
the arbitrary parameter b. 

By performing the integration as in the case of spectrum (12), instead of solution (13) 
in the present case one eventually gets the solution 

= 2afr/37 2 exp [afry 2 ] exp [-am + ] 
' ' 1 — exp [— a(uo + — to>_)] 
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x exp 



c 



sine 



7 2 2 p 2 - (av + 



(21) 



After Mckirmon's, this eq.(21) appears to be the simplest closed-form solution, since 
both of them do not need any recourse to series expansions. In a sense, our solution (21) 
might be regarded as the subluminal analogous of the (superluminal) X-wave solution; 
a difference being that the standard X-shaped solution has a spectrum starting with 0, 
where it assumes its maximum value, while in the present case the spectrum starts at u;_ 
and gets increasing afterwards till u + . More important is to observe that the gaussian 
spectrum has a priori two advantages w.r.t. eq.(20): It may be more easily centered 
around any value ujq of u, and, when increasing its concentration in the surrounding of 
ujq, the spot transverse width does not increase indefinitely, but tends to the spot width 
of a Bessel beam with u = ujq and k z = (ujq — b)/V, at variance with what happens for 
spectrum (20); however, solution (21) is noticeable, since it is really the simplest one. 

Figure [3] show the intensity of the real part of the subluminal pulse corresponding to 
this spectrum, with v = 0.99 c, b = 3 x 10 13 Hz (which result in cu_ = 1.5 x 10 13 Hz and 
u_ = 3 x 10 15 Hz), Au/u + = 1/100 (i.e., a = 100). This is an optical pulse of 0.2 ps. 



( Re T f 



0.9- 

0,8-*,..-' 
0.7-...- 




Figure 2: The intensity of the real part of the subluminal pulse corresponding to spectrum 
(20), with v = 0.99 c, b = 3 x 10 13 Hz (which result in w_ = 1.5 x 10 13 Hz and u;_ = 
3 x 10 15 Hz), Alu/u + = 1/100 (i.e., a = 100). 
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3 A second method for constructing 
subluminal Localized Pulses 



The previous method appears to be very efficient for finding out analytic subluminal LWs, 
but it looses its validity in the limiting case v — > 0, since for v = it is uo^ = uj + and the 
integral in Eq.(3) degenerates, furnishing a null value. By contrast, we are interested also 
in the v = case, since it corresponds to some of the most interesting, and potentially 
useful, LWs: Namely, to the so-called "Frozen Waves", which are new stationary solutions 
to the wave equations, possessing a static envelope. 

It is possible to solve this problem by starting again from Eq.(l), with constraint (2), 
but going on — this time — to integrals over k z . instead of over u. It is enough to write 
relation (2) in the form 

k z = i(w-6) (2') 
for expressing the exact solutions (1) as 

/•/c z max 

V(p,z,t) = expH&t] / dk z S(k z ) Jo(pk p ) exp [i(k z ] , (22) 

J Lmin 



with 



timn = 



and with 



c 1 + 

A;, max = 

cl-/3 



(23) 



k 2 h h 2 
k P 2 = --\ + 2-(3k z + - 2 , (24) 

where quantity ( is still defined according to Eq.(5), always with v < c. 

One can show that the unique exact solution previously known [3] may be rewritten 
in the form of Eq.(24) with S(k z ) = constant. Then, on following the same procedure 
exploited in our first method (previous Section), one can find out new exact solutions 
corresponding to 

S(k z ) = exp l* 2 ^} , (25) 

where 
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K = /c 2 max — k z min , 



by performing the change of variable [analogous, in its finality, to the one in Eq.(7)] 

k z = - 7 2 (s + /5). (26) 
c 

At the end, the exact subluminal solution corresponding to the spectrum (25) results 
to be 



x exp [imx(5] sine 



"\ 2 p 2 +(-~tH + n7c) , (27) 



C 2 V c 



We can again observe that any spectra S(k z ) can be expanded, in the interval /c 2 min < 
max, in a Fourier series: 

OO 2yj- 

S(k z ) = A n ex P [+in—k z ] , (28) 

n=— oo 

with coefficients given now by 

I rk z max 2n 



A n K 



rK z ma,x 2lT 

/ dk z S(k z ) exp [-in— k z \ (29) 

J fc z min -fv 

quantity K having been defined above. 

At the very end of the whole procedure, the general exact solution representing a 
subluminal LW, for any spectra S(k z ), can be eventually written 



x ^ A n exp [m7r/3] sine 



1 



b l 2l 2 p 2 + ( b _ 7 2 ( + nw \ _ (30) 



whose coefficients are expressed in Eq.(29), and where quantity r] is defined as above, in 
Eq.(ll). 

Interesting examples could be easily worked out, as we did at the end of the previous 
Section. 
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4 Stationary solutions with zero-speed envelopes 
("Frozen Waves") 

Here, we shall refer ourselves to the (second) method, expounded in the previous Section. 
Our solution (30), for the case of envelopes at rest, that is, in the case v = [which implies 
( = z], becomes 



b 00 
^(p, z,t) =2 - exp [— ibt] ^ A n sine 



\ 



b 2 (b V 

-ffi+l-z + nv] , (31) 



with coefficients A n given by Eq.(29) with v — 0, so that its integration limits simplify 
into — b/c and b/c, respectively; thus, one gets 



c 



A n = - J^dk z S{k z ) exp [-iny k 2 ] . (29') 

Equation (31) is a new, exact solution, corresponding to stationary beams with a 
static intensity envelope. Let us observe, however, that even in this case one has energy 
propagation, as it can be easily verified from the power flux S = — V\]/ d^/dt (scalar 
case) or from the Poynting flux S = (E A H) (vectorial case: the condition being that ^ 
be a single component, A z , of the vector potential A).[23j For v = 0, eq.(2) becomes 

u = b = lo , 

so that the particular subluminal waves endowed with null velocity are actually monochro- 
matic beams. 

It may be stressed that the present (second) method does yield exact solutions, with- 
out any need of the paraxial approximation, which, on the contrary, is so often used when 
looking for expressions representing beams, like the gaussian ones. Let us recall that, 
when having recourse to the paraxial approximation, the obtained beam expressions are 
valid only when the envelope sizes (e.g., the beam spot) vary in space much more slowly 
than the beam wavelength. For example, the usual expression for a gaussian beam [5*5] 
holds only when the beam spot Ap is much larger than Ao = u;o/(2vrc) = b/(2iTc): so 
that those beams cannot be very much localized. By contrast, our method overcomes 
such problems, since it provides us, as we have seen above, with exact expressions for 
(well localized) beams with sizes of the order of their wavelength. Notice, moreover, that 
the already-known exact solutions — for instance, the Bessel beams — are nothing but 
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particular cases of our solution (31). 



An example: On choosing (with < g_ < q + < 1) the spectral double-step function 



S(k z 



Uo(q+ - q-) 


the coefficients of Eq.(31) become 



A 



ic 



27rnuo(q + — g_ 



for q_u /c < k z < q + u /c 
elsewhere , 



g— iq+irn ^—iq—Trn 



(32) 



(33) 



The double-step spectrum (32), with regard to the longitudinal wave number, corre- 
sponds to the mean value k z = u>o(q + + g_)/2c and to the width Ak z = uj (q + — q-)/c. 
From these relations, it follows that Ak z /k z = 2(q + — q-)/{q+ + (?-)• 

For values of g_ and q + that do not satisfy the inequality Ak z /k z << 1, the resulting 
solution will be a non-paraxial beam. 

Figures 3 show the exact solution corresponding to u = 1.88 x 10 15 Hz (i.e., Ao = 
1 /^m) and to g_ = 0.3, q + = 0.9, which results to be a beam with a spot-size diameter 
of 0.6 fim, and, moreover, with a rather good longitudinal localization. In the case 
of Eqs.(32,33), about 21 terms (—10 < n < 10) in the sum entering Eq.(31) are quite 
enough for a good evaluation of the series. The beam considered in this example is highly 
non-paraxial (with Ak z /k z = 1 ), and therefore couldn't have been obtained by ordinary 
gaussian beam solutions (which are valid in the paraxial regime only)]*} 



Let us now emphasize that a noticeable property of our present method is that it allows 
a spatial modeling even of monochromatic fields (that correspond to envelopes at rest; so 
that, in the electromagnetic cases, one can speak, e.g., of the modeling of "light-fields at 
rest" ). Such a property — rather interesting, especially for applications |62j — was already 
investigated, under different assumptions, in refs. [391 E2 EH] , where the stationary fields 
with static envelopes were called "Frozen Waves" (FW). Namely, in the quoted references, 
discrete superpositions of Bessel beams were adopted in order to get a predetermined 
longitudinal (on-axis) intensity pattern, inside the desired space interval < z < L. In 
other words, in refs.pHl EH ES], the Frozen Waves have been written in the form 

*We are considering here only scalar wave fields. In the case of non-paraxial optical beams, the vector 
character of the field has to be considered 
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m 2 1 




p(nm) 



(b) 



Figure 3: (a) Orthogonal projection of the three-dimensional intensity pattern of the 
beam (a null-speed subluminal wave) corresponding to spectrum (32); (b) 3D plot of the 
field intensity. The beam considered in this example is highly non-paraxial. 



e ~iuJ t e iQz 



N 
n=-N 



i2nirz/c 



(34) 



with 



B n = y [ L dz F{z) e 
L Jo 



■i2nnz/L 



(35) 



quantity |F(^)| 2 being the desired longitudinal intensity shape, chosen a priori. In 
Eq.(35), it is k Pn 2 = n 2 u 2 /c 2 — k zn 2 , and < Q + 2Ntt/L < ujq/c, were we put 
k zn = Q + 2nn/n. As we see from Eq.(34), the FWs have been represented in the past 
in terms of discrete superpositions of Bessel beams. But, now, the method exploited 
in this paper allows us to go on to dealing with continuous superpositions. In fact, the 
continuous superposition analogous to eq.(34) would write 



-iujQt 



ujo/c 



dk z S(k z ) Jo(pk p ) e 



izkz 



(36) 



which, actually, is nothing but our previous Eq.(22) with v = (and therefore ( 
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That is, eq.(36) does just represent a null-speed subluminal wave. To be clearer, let us 
recall that the FWs were expressed in the past as discrete superposition, mainly because 
it was not known at that time how to treat analytically a continuous superposition like 
eq.(36). Only by following the method presented in this work one can eventually extend 
the FW approach [391 ES EH] to the case of integrals: without numerical simulations, but 
in terms once more of analytic solutions. 

Indeed, the exact solution of Eq.(36) is given by Eq.(31), with coefficients (29'). And 
one can choose the spectral function S(k z ) in such a way that \1/ assumes the on-axis pre- 
chosen static intensity pattern |F(z)| 2 . Namely, the equation to be satisfied by S(k z ), to 
such an aim, comes by associating eq.(36) with the requirement |^(p = 0, z, t)\ 2 = \F(z)\ 2 , 
which entails the integral relation 

r' C dk z S(k z ) e izk * = F(z) . (37) 

J-ljo/c 

Equation (37) would be trivially soluble in the case of an integration between — oo and 
+oo, since it would merely be a Fourier transformation; but obviously this is not the 
case, because its integration limits are finite. Actually, there are functions F(z) for which 
Eq.(37) is not soluble, in the sense that no spectra S{k z ) exist obeying the last equation. 
Namely, if we consider the Fourier expansion 

dk z S(k z ) e izk > , 

-oo 

when S(k z ) does assume non-negligible values outside the interval —uj /c < k z < uj /c, 
then in Eq.(37) no S(k z ) can forward that particular F(z) as a result. 

However, way-outs can be devised, such that one can can nevertheless find out a 
function S{k z ) that approximately (but satisfactorily) complies with our Eq.(37). 

The first way-out consists in writing S(k z ) in the form 

S(k z ) = I £ F (^) e ~i2n^/K (3g) 

where, as before, K = 2ujq/c. Then, one can easily verify Eq.(38) to guarantee that the 
integral in Eq.(37) yields the values of the desired F(z) at the discrete points z = 2nir/K. 
Indeed, the Fourier expansion (38) is already of the same type as Eq.(28), so that in this 
case the coefficients A n of our solution (31), appearing in Eq.(29'), do simply become 
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(39) 



This is a powerful way for obtaining a desired longitudinal (on-axis) intensity pattern, 
especially for very small spatial regions, because it is not necessary to solve any integral 
to find out the coefficients A n , which by contrast are given directly by Eq.(39). 

Figures [3] depict some interesting applications of this method. A few desired longitu- 
dinal intensity patterns |F(z)| 2 have been chosen, and the corresponding Frozen Waves 
calculated by using Eq.(31) with the coefficients A n given in Eq.(39). The desired patterns 
are enforced to exist within very small spatial intervals only, in order to show the capability 
of our method to model the field intensity shape also under such strict requirements. 

In the four examples below, we considered a wavelength A = 0.6 /zm, which corre- 
sponds to uj = b = 3.14 x 10 15 Hz. 

The first longitudinal (on-axis) pattern considered by us is that given by 



( e a{z ~ z) for < z < Z 

F{z) = 

[ elsewhere 

i.e., a pattern with an exponential increase, starting from z = till z = Z. The chosen 
values of a and Z are Z = 10 /im and a = 3/Z. The intensity of the corresponding Frozen 
Wave is shown in Fig. (4a). 

The second longitudinal pattern (on-axis) taken into consideration is the gaussian 
one, given by 

f e" 9( t )2 for -Z<z<Z 

F{z) = 

{ elsewhere 

with q = 2 and Z = 1.6 /im. The intensity of the corresponding Frozen Wave is shown in 
Fig. (4b). 

In the third example, the desired longitudinal pattern is supposed to be a super- 
gaussian: 
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F(z) 



e -«(f ) 2m for - Z < z < Z 



elsewhere 



where m controls the edge sharpness. Here we have chosen q = 2, m = 4 and Z = 2 //m. 
The intensity of the Frozen Wave obtained in this case is shown in Fig. (4c). 

Finally, in the fourth example, let us choose the longitudinal pattern as being the 
zero-order Bessel function 



!Jo(q z) for — Z < z < Z 
elsewhere 

with q = 1.6 x 10 6 m _1 and Z = 15 /im. The intensity of the corresponding Frozen Wave 
is shown in Fig.(4d). 

Let us observe that, of course, any static envelopes of this type can be easily trans- 
formed into propagating pulses by the mere application of Lorentz transformations. 



Another way-out exists for evaluating S(k z ), based on the assumption that 

S(k z ) ~ S(k z ) , (40) 

which consitutes a good approximation whenever S(k z ) assumes negligible values outside 
the interval [—uq/c, ujq/c]. In such a case, one can have recourse to the method associated 
with Eq.(28) and expand S(k z ) itself in a Fourier series, getting eventually the relevant 
coefficients A n by Eq.(29). Let us recall that it is still K = fc z max — fc 2 min = 2uj /c. 

It may be interesting to call attention to the circumstance that, when constructing 
FWs in terms of a sum of discrete superpositions of Bessel beams (as it was done by us in 
refs. [39| I4TR l59| 162]). it was easy to obtain extended envelopes like, e.g., "cigars": where 
easy means having recourse to few terms of the sum. By contrast, when we construct 
FWs — following this Section — as continuous superpositions, then it is easy to get higly 
localized (concentrated) envelopes. Let us explicitly mention, moreover, that the method 
presented in this Section furnishes FWs that are no longer periodic along the z-axis (a 
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(a) (b) 




Figure 4: Frozen Waves with the on-axis longitudinal field pattern chosen as: (a) Expo- 
nential; (b) Gaussian; (c) Super-gaussian; (d) Zero order Bessel function 

situation that, with our old method, [39l HDJ [59] was obtainable only when the periodicity 
interval tended to infinity). 



5 The role of Special Relativity, and of Lorentz 
Transformations 

Strict connections exist between, on one hand, the principles and structure of Special 
Relativity and, on the other hand, the whole subject of subluminal, luminal, superluminal 
Localized Waves, in the sense that it is expected since long time that a priori they are 
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transformable one into the other via suitable Lorentz transformations (cf. refs.[T| IW[ IM], 
besides work of ours in progress)^ 

Let us first confine ourseves to the cases faced in this paper. Our subluminal localized 
pulses, that may be called "wave bullets", behave as particles: Indeed, our subluminal 
pulses [as well as the luminal and superluminal (X-shaped) ones, that have been amply 
investigated in the past literature] do exist as solutions of any wave equations, ranging 
from electromagnetism and acoustics or geophysics, to elementary particle physics (and 
even, as we discovered recently, to gravitation physics). From the kinematical point of 
view, the velocity composition relativistic law holds also for them. The same is true, 
more in general, for any localized waves (pulses or beams). 

Let us start for simplicity by considering, in an initial reference-frame O, just a (y- 
order) Bessel beam: 

0, z, t) = J v {pk p ) e** e izk * e-™ . (41) 

In a second reference- frame O', moving with respect to (w.r.t.) O with speed u — along the 
positive z-axis and in the positive direction, for simplicity's sake — , it will be observed [dD] 
the new Bessel beam 

V(j/,<ff,z?,lf) = J v {p'k',s) e** e iz ' k '*' e'™' 1 ' , (42) 

obtained by applying the appropriate Lorentz transformation (a Lorentz "boost") with 
7 = [y/l-ut/c?]- 1 : 

k'p' = k p ; k' z i = "i{k z — uu/c 2 ); uj' = j(uj — uk z ) ; (43) 
this can be easily seen, e.g., by putting 

p = p; z = 7(2' + ut'); t = 7(t' + uz'/c 2 ) (44) 

directly into Eq.(42). 

Let us now pass[60j to subluminal pulses. We can investigate the action of a Lorentz 
transformation (LT), by expressing them either via the first method (Section 2) or via 

§After the completion of this work, we realized that we had overlooked an important and interesting 
paper by Saari et al.|60j. wherein the relativistic connections between the LWs are already, deeply and 
in general, investigated in terms of suitable LTs: we therefore inserted in this article the necessary 
quotations. We are actually glad in calling attention to ref. [60], also because the inspiring philosophy 
— which in part goes back to papers like refs.[U[TH] — is shared by us too. 
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the second one (Section 3). Let us consider for instance, in the frame O, a v -speed 
(subluminal) pulse, given by Eq.(3) of our Section 2. When we go on to a second observer 
O' moving with the same speed v w.r.t. frame O, and, still for the sake of simplicity, passing 
through the origin O of the initial frame at time t = 0, the new observer O' will see the 
pulse [60J 

*(f/,z',1?) = e- it,UJ '° r + duo S(u) Mp'k'p,) e iz ' k '*' , (45) 

Jul- 

with 

k' z i = , j~ 1 oj/v — 'jb/v; uJ = 76; k' p> = oj'o/c 2 — k' z / 2 , (46) 

as one gets from the Lorentz transformation in Eq.(43), or in Eq.(44), with u = v [and 
7 given by Eqs.(9)]. Notice that k' ' z > is a function of u, as expressed by the first one of 
the three relations in the previous Eqs.(46); and that here to' is a constant. 

If we explicitly insert into Eq.(45) the relation uj = j(v k' z > + 76), which is nothing 
but a re-writing of the first one of Eqs.(46), then Eq.(45) becomes [60] 

*(p', z', t') = 1V e~ u ^ / dk' z , S(k' z ,) Mp'k'p,) e lz k *> , (47) 
where S is expressed in terms of the previous function S(u), entering Eq.(45), as follows: 

S(k' zl ) = S( 1V k' z , +1 2 b) . (48) 

Equation (47) describes monochromatic beams with axial symmetry (and does coincide 
also with what derived within our second method, in Section 3, when posing v — 0). 

The remarkable conclusion is that a subluminal pulse, given by our Eq.(3), which 
appears as a v -speed pulse in a frame O, will appear [60J in another frame O' (traveling 
w.r.t. observer O with the same speed v in the same direction z) just as the monochromatic 
beam in Eq. (47) endowed with angular frequency u' = 76, whatever be the pulse spectral 
function in the initial frame O: even if the kind of monochromatic beam one arrives to 
does of course depencj^on the chosen S(u). The vice- versa is also true, in general. 

Let us set forth explicitly an observation that hasn't been noticed in the existing 
literature yet. Namely, let us mention that, when starting not from Eq.(3) but from 

^One gets in particular a Bessel-type beam when S is a Dirac's delta-function: S(oj) = 5(uj — uq)- 
Moreover, let us notice that, on applying a LT to a Bessel beam, one obtains another Bessel beam, with 
a different axicon-angle. 
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the most general solutions which — as we have already seen — are sums of solutions (3) 
over the various values b m of b, then a Lorentz transformation will lead us to a sum of 
monochromatic beams: actually, of harmonics (rather than to a single monochromatic 
beam). In particular, if one wants to obtain a sum of harmonic beams, one has to apply 
a LT to more general subluminal pulses. 




Figure 5: Let us consider an object that is intrinsically spherical, i.e., that is a sphere in 
its rest-frame (Panel (a)). After a generic subluminal LT along x, i.e., under a subluminal 
x-boost, it is predicted by Special Relativity (SR) to appear as ellipsoidal due to Lorentz 
contraction (Panel (b)). After a superluminal x-boost [T6l 1651 164"] (namely, when this 
object moves|21J with superluminal speed V), it is predicted by SR — in its non-restricted 
version (ER) — to appear [TJ as in Panel (d), i.e., as occupying the cylindrically symmetric 
region bounded by a two-sheeted rotation hyperboloid and an indefinite double cone. The 
whole structure, according to ER, is expected to move rigidly and, of course, with the 
speed V, the cotangent square of the cone semi-angle being (V/c) 2 — 1. Panel (c) refers 
to the limiting case when the boost-speed tends to c, either from the left or from the right 
(for simplicity, a space axis is skipped). It is remarkable that the shape of the localized 
(subluminal and superluminal) pulses, solutions to the wave equations, appears to follow 
the same behavior; this can have a role for a better comprehension even of de Broglie and 
Schroedinger wave- mechanics. The present figure is taken from refs.[T| [To]. See also the 
next Figure. 

Let us add that also the various superluminal localized pulses get transformed jfiT)] one 
into the other by the mere application of ordinary LTs; while it may be expected that the 
subluminal and the superluminal LWs are to be linked (apart from some known technical 
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difficulties, that require a particular caution) by the superluminal Lorentz "transforma- 
tions" expounded long ago, e.g., in refs.pilj ESI EH H] and refs. therein Let us recall 
once more that, in the years 1980-82, special relativity, in its non-restricted version, pre- 
dicted that, while the simplest subluminal object is obviously a sphere (or, in the limit, a 
space point), the simplest superluminal object is on the contrary an X-shaped pulse (or, 
in the limit, a double cone): cf. Figj5j taken from refs.[U [16]. The circumstance that 
the pattern of the localized solutions to the wave equations does meet this prediction is 
rather interesting, and is expected to be useful — in the case, e.g., of elementary particles 
and quantum physics — for a deeper comprehension of de Broglie's and Schroedinger's 
wave mechanics. With regard to the fact that the simplest subluminal LWs, solutions to 
the wave equation, are "ball-like", let us depict by Figs|6j in the ordinary 3D space, the 
general shape of the Mackinnon's solutions as expressed by Eq.(10) for v << c: In such 
figures we graphically represent the field iso-intensity surfaces, which in the considered 
case result to be (as expected) just spherical. 

We have also seen, among the others, that, even if our first method (Section 2) 
cannot yield directly zero-speed envelopes, such envelopes "at rest", in Eq.(31), can be 
however obtained by applying a n-speed LT to eq.(16). In this way, one starts from 
many frequencies [Eq.(16)] and ends up with one frequency only [Eq.(31)], since b gets 
transformed into the frequency of the monochromatic beam. 



6 Non-axially symmetric solutions: The case of 
higher-order Bessel beams 

Let us stress that till now we have paid attention to exact solutions representing 
axially-symmetric (subluminal) pulses only: that is to say, to pulses obtained by suitable 
superpositions of zero-order Bessel beams. 

It is however interesting to look also for analytic solutions representing non-axially 
symmetric subluminal pulses, which can be constructed in terms of superpositions of v- 

I'One should pay attention that, as we were saying, the topic of superluminal LTs is a delicate [TBI 
1551 ?, UJ one, at the extent that the majority of the recent attempts to re-address this question and its 
applications seem to be defective (sometimes they do not even keep the necessary covariance of the wave 
equation itself). 
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Figure 6: In the previous Figure we have seen how SR, in its non-restricted version (ER), 
predicted [U HE] that, while the simplest subluminal object is obviously a shere (or, in the 
limit, a space point), the simplest superluminal object is on the contrary an X-shaped 
pulse (or, in the limit, a double cone). The circumstance that the Localized Solutions to 
the wave equations do follow the same pattern is rather interesting, and is expected to 
be useful — in the case, e.g., of elementary particles and quantum physics — for a deeper 
comprehension of de Broglie's and Schroedinger's wave mechanics. With regard to the 
fact that the simplest subluminal LWs, solutions to the wave equations, are "ball-like" , let 
us depict by these figures, in the ordinary 3D space, the general shape of the Mackinnon's 
solutions as expressed by Eq.(lO), numerically evaluated for v « c. In figures (a) and 
(b) we graphically represent the field iso-intensity surfaces, which in the considered case 
result to be (as expected) just spherical. 

order Bessel beams, with v a positive integer {y > 0). This can be attempted both in 
the case of Sect. 2 (first method), and in the case of Sect. 3 (second method). 

For brevity's sake, let us take only the first methos (Sect. 2) into consideration. One is 
immediately confronted with the difficulty that no exact solution is known for the integral 
in Eq.(8) when J (.) is replaced with J v (.). 

One can overcome this difficulty by following a simple method, which will allow us to 
obtain "higher-order" subluminal waves in terms of the axially-symmetric ones. Indeed, 
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it is well-known that, if ^(x, y, z, t) is an exact solution to the ordinary wave equation, 
then d n ^> /dx n and d n ^/dy n are also exact solutions]^] One should notice that, on 
the contrary, when working in cylindrical co-ordinates, if \I/(p, 0, z,t) is a solution to the 
wave equation, /dp and <9\I//<90 are not solutions, in general. Nevertheless, it is not 
difficult at all to reach the noticeable conclusion that, once \I/(p, 0, z,t) is a solution, 
then also 

is an exact solution. For example, for an axially-symmetric solution 



of the type \I/ = Jo(k p p) exp[ik z ] exp[—iut], equation (49) yields \I/ = 
—k p Jx(kpp) exp[i0] exp[ik z ] exp[—iut], which is actually another analytic solution. 

In other words, it is enough to start for simplicity from a zero-order Bessel beam, 



and to apply Eq.(49), successively, v times, in order to get as a new solution \I/ = 
(—k p ) u J u (k p p) exp[iu(j)} exp[ik z ] exp[— iut], which is a i/-order Bessel beam. 

In such a way, when applying v times Eq.d49j) to the (axially-symmetric) subluminal 
solution *&(p, z, t) in Eqs.(16,15,14) [obtained from Eq.(3) with spectral function S(u)], 
we are able to get the subluminal non-axially symmetric pulses ^> u (p,4>, z,t) as new 
analytic solutions, consisting as expected in superpositions of //-order Bessel beams: 

tf„(p, 0, z, t) = P du S'(uj) J v {k p p) e 1 ^ e lk * z e"*" 1 , (50) 



with kp{uS) given by Eq.(4), and quantities S'(uS) = (—k p (uj)) u S(uj) being the spectra 
of the new pulses. If S(u) is centered at a certain carrier frequency (it is a gaussian 
spectrum, for instance), then S'(u) too will approximately be of the same type. 

Now, if we wish the new solution \I/„(p, 0, z, t) to possess a pre-defined spectrum 
S'(uj) = F(u), we can first take Eq.(3) and put S(u) = F(u)/(—k p (u)) u in its solution 
(16), and afterwards apply to it, v times, the operator U = exp[i(p] [d/dp + (i/p)d/d(p)} : 
As a result, we will obtain the desired pulse, ^ u (p, 0, z, t), endowed with S'(u) = F(u). 

An example: 

On starting from the subluminal axially-symmetric pulse ^(p, z, t), given by Eq.(16) with 
the gaussian spectrum (17), we can get the subluminal, non-axially symmetric, exact 



**Let us mention that even d n ^/dz n and d ni l> /dt n will be exact solutions. 
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solution x I / i(p, <p, z, t) by simply calculating 



dp 



,i<t> 



(51) 



which actually yields the "first-order" pulse \l/i(p, <p, z, t), which can be more compactly 
written in the form: 



0,^,0 = 2-^7 2 exp 
c 



i- /?7 2 ?7 
c 



7T 

^ A„ exp [in-] Vin 

i=-oo P 



with 



where 



ip ln (p,(f),r),C) = ^7 2 P ^ 3 [Z cos Z - sin Z] e* 



(52) 



(53) 



Z = 



\ 



7 2 p 2 + - 7 2 ( + wr . 



(54) 



This exact solution, let us repeat, corresponds to superposition (50), with S'(uj) = 
k p (uj)S(uj), quantity S(uS) being given by Eq.(17). It is represented in Figure [7j The 
pulse intensity has a "donut-like" shape. 




C (nm) 



Figure 7: Orthogonal projection of the field intensity corresponding to the higher order 
subluminal pulse represented by the exact solution Eq.(51), quantity \l/ being given by 
Eq.(16) with the gaussian spectrum (17). The pulse intensity happens to have this time 
a "donut"-like shape. 



28 



7 Conclusions 



Like in the well-known superluminal case, the (more orthodox, in a sense) subluminal 
Localized Waves can be obtained by suitable superpositions of Bessel beams. They had 
been till now almost neglected, however, for the mathematical difficulties met in finding 
out analytic expressions for them, difficulties associated with the fact that the superposi- 
tion integral runs over a finite interval. We have re-addressed in this paper the question of 
such subluminal LWs, obtaining by contrast, and in a simple way, non-diffracting sublu- 
minal pulses as exact analytic solutions to the wave equations. Such new ideal subluminal 
solutions — which propagate without distortion in any homogeneous linear medium — have 
been herein obtained for arbitrarily chosen frequencies and bandwidths, avoiding in par- 
ticular any recourse to the non-causal components so frequently plaguing the previously 
known localized waves. 

Indeed, only one closed-form subluminal LW solution, ip cil to the wave equations was 
known|3j: the one obtained by choosing in the relevant integration a constant weight- 
function S(u); while all other solutions had been got in the past only via numerical 
simulations. By contrast, we have shown above, for instance, that a subluminal LW can 
be obtained in closed form by adopting any spectra S(u) that be expansions in terms of 
ip c f. In fact, the initial disadvantage of the subluminal case, that of being associated with 
a a limited bandwidth, has been turned into an advantage, since in the case of "truncated" 
integrals — at variance with the superluminal case — the spectrum S(u) can be expanded 
in a Fourier series. 

More in general, it has been shown in what precedes how can one arrive at exact 
solutions both by integration over the Bessel beams' angular frequency lu, and by in- 
tegration over their longitudinal wavenumber k z . Both methods are expounded in this 
paper. The first one appears to be powerful enough; we have studied the second method 
as well, however, since it allows dealing even with the limiting case of zero-speed solutions 
(thus furnishing a new way, in terms of continuous spectra, for obtaining the so-called 
"Frozen Waves", so promising also from the point of view of applications). We have 
briefly treated the case, moreover, of non axially-symmetric solutions, namely, of higher 
order Bessel beams. 

At last, attention has been always paid to the role of Special Relativity, and to the 
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fact that the localized waves are transformed one into the other by suitable Lorentz 
Transformations. Moreover, our results seem to show that in the subluminal case the 
simplest LW solutions are (for v « c) "ball"-like, as expected since longjl] on the mere 
basis of special relativity [16]. Actually, already in the years 1980-82 it had been predicted 
that, if the simplest subluminal object is a sphere (or, in the limit, a space point), then 
the simplest superluminal object is an X-shaped pulse (or, in the limit, a double-cone); 
and viceversa: Cf. Fig.5, taken from ref. pp. It is rather interesting that the same pattern 
appears to be followed by the localized solutions of the wave equations. For the subluminal 
case, see. e.g., Figs. 6. 

The localized pulses, endowed with a finite energy, or merely truncated, will be con- 
structed in another paper. 

In the present work we have fixed our attention especially on electromagnetism and 
optics: but results of the same kind are valid whenever an essential role is played by a 
wave-equation (like in acoustics, seismology, geophysics, gravitation, elementary particle 
physics, etc.). 
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